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Brief  Outline  of  Research  Findings 


Section  II  of  this  report  indicates  the  publication  status  of  the  research  work  of  mem¬ 
bers  of  the  Applied  Mathematics  Group.  A  number  of  papers  previously  accepted  have 
been  published,  others  previously  submitted  have  been  accepted,  and  a  number  of  new 
ones  have  been  submitted. 

Among  the  latter  are  two  papers  by  Professor  Keller  and  his  former  student  Dr.  Kevin 
C.  Nunan.  They  calculated  the  effective^  elastic  behaviour  of  a  composite  of  rigid  spheres 
located  at  lattice  points  in  an  elastic  matrix.  They  also  calculated  the  effective  viscosity  of 
a  similar  arrangement  of  rigid  spheres  in  a  viscous  fluid.  The  novelty  of  this  work  is  that 
it  covers  all  concentrations  from  zero  to  close  packing. 

■  Ar 

Professors  Geer  and  Keller  completed  their  work  on  the  eigenvalues  of  slender  cavities, 

which  provides  a  simple  method  for  calculating  such  eigenvalues  accurately.  Their  results 
agreed  very  well  with  previous  numerical  calculations  for  prolate  spheroids. 

Dr.  Margaret  Cheney  has  completed  two  new  works  on  inverse  scattering  in  two  di¬ 
mensions.  This  provides  a  further  step  towards  making  inverse  scattering  theory  applicable 
to  real  physical  problems. 
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II.  PUBLICATIONS  OF  APPLIED  MATHEMATICS  GROUP 


50.  J.  B.  Keller 


Optimum  inspection  policies 

Pub:  Management  Sci.,  Z8,  447-450,  1982. 


51.  P.  S.  Hagan  Travelling  and  stacked  travelling  wave 

solutions  of  parabolic  equations  J 

Pub:  SIAM  J.  Math.  Anal.,  13,  717-738,  1982. 

52.  J.-M.  Vanden-Broeck  Contact  problems  involving  the  flow  past  an 

inflated  aerofoil 

Pub;  J.  Appl.  Mech.,  49,  263-265,  1982. 


53.  J.-M.  Vanden-Broeck  Nonlinar  two-dimensional  sail  theory 

Pub:  Phys.  Fluids,  25,  420-423,  1982. 

55.  P.  S.  Hagan  Spiral  waves  in  reaction  diffusion  equations 

Pub:  SIAM  J.  Appl.  Math.,  42,  762-786,  1982. 


56.  J.-M.  Vanden-Broeck  Parabolic  approximations  for  ship  waves  and 

J.  B.  Keller  wave  resistance 

Pub:  Proceedings  of  the  Third  International 

Conference  on  Numerical  Ship  Hydrodynamics, 
Paris,  France,  June  16-19,  1981. 


57.  A.  Jeffrey 
T.  Kawahara 


Asymptotic  Methods  in  Nonlinear  Wave  Problems 
Pub:  Pitman  Publishing,  Ltd.,  London,  1982. 


58.  M.  J.  Miksis 

J.-M.  Vanden-Broeck 
J.  B.  Keller 


Rising  Bubbles 

Pub:  J.  Fluid  Mech.,  123,  31-42,  1982. 


59.  J.  C.  Neu 


60.  J.  B.  Keller 
M.  J.  Miksis 


61.  A.  Jepson 
A.  Spence 


Resonantly  interacting  waves 

Pub:  SIAM  J.  Appl.  Math. ,  43,  I,  141-156,  1983.' 
Surface  tension  driven  flows 

Pub:  SIAM  J.  Appl.  Math.,  43,  II,  268-277,  1983. 

Folds  in  solutions  of  two  parameter  systems 
and  their  calculation;  Part  I 

Pub :  Stanford  Univ.  Numer.  Anal.  Report, 

82-02. 


62.  J.  B.  Keller 


Time-dependent  queues 

Pub:  SIAM  Rev.,  24,  401-412,  1982 


63.  R.  E.  Caflisch 
B.  Nicolaenko 


64.  P.  S.  Hagan 

R.  E.  Caflisch 

J.  B.  Keller 


65.  R.  E.  Caflisch 


66. 


J.  B.  Keller 
J.-M.  Vanden-Broeck 


67.  R.  E.  Caflisch 


68.*  M.  S.  Falkovitz 
M  Seul 
H.  L.  Frisch 
H.  M.  McConnell 


69.  R.  E.  Caflisch 


70.  P.  F.  Rhodes-Robinson 


71.  P.  F.  Rhodes-Robinson 


Shock  profile  solutions  of  the  Boltzmann  equation 
Pub :  Comm.  Math.  Phys.,  86,  161-194,  1982. 

Arrow' s  model  of  optimal  pricing,  use  and 
exploration  of  under tain  natural  resources 

Sub:  Econometrics 

Radiation  transport  in  a  hot  plasma 
Acc:  SIAM  J.  Appl.  Math.,  in  press. 

Jets  rising  and  falling  under  gravity 
Pub:  J.  Fluid  Mech.,  124,  335-345,  1982. 

Fluid  dynamics  and  the  Boltzmann  equation 
Pub:  Stud.  Stat.  Mech.,_5,  194-223,  1983. 

Theory  of  periodic  structures  in  lipid 
bilayer  membranes 

Pub:  Proc.  Nat.  Acad.  Sci.,  7_9,  3918,  1982. 

The  fluid-dynamic  limit  of  a  model  Boltzmann 
equation  in  the  presence  of  a  shock 

Pub:  Institute  National  de  Recherche  en 

Informatique  et  en  Automatique,  INRIA 
No.  81,  June  1981,  1-34. 


On  the  short  surface  waves  due  to  half-immersed 
circular  cylinder  oscillating  on  water  of 
infinite  depth 

Pub:  Proc.  Royal  Soc.  London  A,  384,  333-357,  1982. 

Note  on  the  reflexion  of  water  waves  at  a 
wall  in  the  presence  of  surface  tension 

Pub:  Math.  Proc.  Cambridge  Philosophical  Soc., 

92,  369-373,  1982. 
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72.  P.  F.  Rhode s-Robinson 


R.  E.  Caflisch 
G.  C.  Papanicolaou 


R.  E.  Caflisch 
G.  C.  Papanicolaou 


75.  R.  E.  Caflisch 
B.  Nicolaenko 


76.  J.  H.  Maddocks 


M.  S.  Falkovitz 
L.  A.  Segel 


M.  S.  Falkovitz 
J.  L.  Frisch 


79.  M.  S.  Falkovitz 

H.  L.  Frisch 

J.  B.  Keller 


80.*  H.  S.  Falkovitz 
J.  B.  Keller 


81.  A.  Jeffrey 
J.  Mvungi 


*Not  supported  by  AFOSR  or  OMR. 


On  the  generation  of  water  waves  at  an 
inertial  surface 

Acc:  J.  Australian  Math.  Soc.  B,  in  press. 


Dynamic  theory  of  suspensions  with  Brownian 
effects 

Acc:  SIAM  J.  Appl.  Math.,  in  press. 


Instability  in  settling  of  suspensions  due 
to  Brownian  effects 

Pub :  Proceedings  of  conference  Two-Phase  Flow. 


Shock  waves  and  the  Boltzmann  equation 

Pub:  Proceedings  NSF-AMS  conference  non-linear 
PDE. 


Restricted  quadratic  forms  and  their  application 
to  bifurcation  and  stability  in  constrained 
variational  principles 

Sub:  SIAM  J.  Appl.  Math. 


Spatially  inhomogeneous  polymerization  in 
unstirred  bulk 

Pub:  SIAM  J.  Appl.  Math.,  42,  386-416,  1983. 


The  scale  of  non-homogeneity  as  defined  by 
diffusion  measurements 

Pub:  Journal  of  Membrane  Science,  10,  61,  1982. 


Optimal  catalyst  distribution  in  a  membrane 
Acc;  Chem.  Eng.  Sci.,  in  press. 

Crawling  of  Worms 

Acc:  J.  Theor.  Biol.,  in  press. 

The  random  choice  method  and  the  free-surface 
water  profile  after  the  collapse  of  a  dam  wall 

Pub:  Wave  Motion,  4,  381-389,  1982. 
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82. 


J.  G.  Watson 
E.  L.  Reiss 


A  statistical  theory  for  imperfect  bifurcation 
Pub:  SIAM  J.  Appl .  Math.,  _42,  135-148,  1982. 


83.  J.  G.  Watson 
J.  B.  Keller 

84.  M.  I.  Weinstein 


Reflection,  scattering  and  absorption  of  acoustic 
waves  by  rough  surfaces 

Acc  :  J*  Acoust.  Soc.  Am. 

Global  existence  for  a  generalized  Korteweg  - 
de  Vries  equation 

Sub:  SIAM  J.  Math.  Anal. 


85.  M.  I.  Weinstein 


86.  M.  Cheney 


Nonlinear  Schrodinger  equations  and  sharp 

interpolation  estimates 

Pub;  Commun.  Math.  Phys.,  87,  567-576,  1983. 


Two-dimensional  scattering:  the  number  of 
bound  states  from  scattering  data 

Sub:  J.  Math.  Phys. 


87.  L.  L.  Bonilla 
A.  Linan 


Relaxation  oscillations,  pulses,  and  travelling 
waves  in  the  diffusive  Volterra  delay-differential 
equation 

Acc:  SIAM  J.  Appl.  Math.,  in  press. 


88.  P.  F.  Rhodes-Robinson  Note  on  the  effect  of  surface  tension  on  water 

waves  at  an  inertial  surface 

Pub:  J.  Fluid  Mech.,  125,  375-377,  1982. 


89.  J.  B.  Keller 
J.  Hunter 


Weak  shock  diffraction 
Acc:  Wave  Motion,  in  press. 


90.  J.  B.  Keller 
J.  Hunter 


Weakly  nonlinear  high  frequency  waves 

Pub:  Comm.  Pure  Appl.  Math.,  36,  547-569,  1983. 


91.  J.  B.  Keller 
J.  C.  Neu 


Asymptotic  analysis  of  a  viscous  Cochlear  model 
Sub :  J.  Acoust.  Soc.  Am. 


92.  J.  B.  Keller 

J.-M.  Vanden-Broeck 


Parabolic  approximations  for  ship  waves  and  wave 
resistance 

Pub:  Proc.  3rd  Inti.  Conf.  on  Numerical  Ship 

Hydrodynamics, Paris,  France,  June  16-19,1981, 


93.  A.  Spence 
A.  Jepson 


J.  B.  Keller 
R.  Burridge 


J.  B.  Keller 


J.  B.  Keller 
A.  S.  Whittemore 


J.  B.  Keller 
J.  F.  Geer 


J.  B.  Keller 
R.  Voronka 


The  numerical  computation  of  turning  points  of 
nonlinear  equations 

Pub:  Treatment  of  Integral  Equations  by  Numerical 
Methods,  169-183,  London,  1982. 


Biot’s  poroelasticity  equations  by  homogenization 

Pub:  Springer  Lecture  Notes,  in  Macroscopic 

Properties  of  Disordered  Media,  NY,  51-57,1982. 


Capillary  waves  on  a  vertical  jet 
Acc:  J.  Fluid  Mech.,  in  press. 

Survival  estimation  using  censored  data 
Sub:  J.  Royal  Statist.  Soc.,  Series  B 

Eigenvalues  of  slender  cavities  and  waves  in 
slender  tubes 

Acc:  J.  Acoust.  Soc.  An.,  in  press. 

Valuation  of  stocks  and  options 
To  be  submitted. 


M.  Cheney 


Inverse  scattering  in  dimension  two 
Acc :  J.  Math.  Phys.,  in  press. 


K.  C.  Nunan 
J.  B.  Keller 


Effective  viscosity  of  a  periodic  suspension 
Sub:  J.  Fluid  Mech. 


K.  C.  Nunan 
J.  B.  Keller 


Effective  elasticity  Tensor  of  a  Periodic  Composite 
Acc:  J.  Mech.  Phys.  Solids,  in  press. 


102.  J.  B.  Keller 


Breaking  of  liquid  films  and  threads 
Acc:  Phys.  Fluids,  in  press. 


103.  J.  B.  Keller 
G.  R.  Verma 


Hanging  rope  of  minimun  elongation 
Acc:  SIAM  Rev. 
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M.  Cheney 
S.  Coen 
A.  Weglein 


J.  B.  Keller 


106.  S.  Venakides 


107.  J.  H.  Maddocks 


108.*  J.  B.  Keller 


109.  J.  B.  Keller 

M.  S.  Falkovitz 


Velocity  &  density  of  a  two-diraensional  acoustic 
medium  from  point  source  surface  data 

Sub:  Phys.  Rev.  Lett. 


Probability  of  a  shutout  in  racquetball 
Acc  i  SIAM  Rev. 

The  Zero  Dispersion  Limit  of  the  Korteweg-de  Vries 
Equation  for  Initial  Potentials  with  Non-trivial 
Reflection  Coefficient 

Sub :  Comm.  Pure  Appl.  Math. 

Stability  of  Nonlinearly  Elastic  Rods 
Acc:  Arch.  Rat.  Mech.  Anal.,  in  press. 

Genetic  Variability  Due  to  Geographical  Inhomogeneity 

Sub:  J.  Math.  Biol. 

Precipitation  pattern  formation 
In  preparation. 
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III.  ABSTRACTS  OF  MANUSCRIPTS  SUBMITTED  OVER  REPORT  PERIOD 


1.  WEAK  SHOCK  DIFFRACTION,  by  J.  Hunter  and  J.B.  Keller. 

The  diffraction  of  a  weak  shock  by  a  rigid  wedge  is  analyzed  theoretically  via  the  theory 
of  weakly  nonlinear  geometrical  acoustics,  which  is  the  same  as  Whitham’s  nonlinearization 
technique.  First  linear  and  weakly  nonlinear  geometrical  acoustics  are  explained.  Then 
the  linear  acoustics  results  for  weak  shock  diffraction  by  a  wedge  are  presented.  Next  these 
results  are  modified  according  to  the  principles  of  weakly  nonlinear  geometrical  acoustics. 
The  results  show  that  the  compressive  diffracted  wavefronts  of  linear  acoustics  are  actually 
shocks,  and  their  positions  and  strengths  are  found.  The  infinite  gradients  of  the  linear 
acoustics  rarefaction  waves  are  found  to  be  finite  but  discontinuous  gradients.  Finally 
the  results  arc  specialized  to  a  shock  hitting  a  right-angled  wedge,  a  shock  corning  off  a 
right-angled  wedge,  and  a  shock  hitting  a  thin  semi-infinite  screen. 


2.  WEAKLY  NONLINEAR  HIGH  FREQUENCY  WAVES,  by  J.  Hunter  and  J.B.  Keller. 

We  shall  present  a  method  for  finding  small-amplitude,  high-frequency  wave  solutions 
of  hyperbolic  systems  of  quasilinear  partial  differential  equations.  The  method  applies  in 
any  number  of  space  dimensions,  and  reduces  to  “geometrical  optics”  when  the  equations 
are  linear.  Therefore  we  call  it  weakly  nonlinear  geometrical  optics. 

The  method  yields  the  first  terms  in  an  asymptotic  expansion  of  a  solution  with 
respect  to  a  small  parameter  £.  This  parameter  is  a  measure  of  both  the  amplitude  and 
the  period  of  a  wave.  The  solution  consists  of  a  slowly  varying  mean  and  a  superposition 
of  small-amplitude,  high-frequency  waves.  Each  of  these  waves  undergoes  distortion  due 
to  nonlinear  selfinteraction.  When  the  distortion  leads  to  multi-valuedness,  shocks  are 
introduced  to  maintain  single-valuedness.  For  conservative  systems,  the  shock  locations 
are  found  to  obey  an  equal  area  rule.  The  interaction  among  the  various  waves  is  a  higher- 
order  effect  which  is  negligible  unless  a  certain  resonance  condition  bolds.  The  present 
method  does  not  apply  in  the  resonant  case. 


3.  CAPILLARY  WAVES  ON  A  VERTICAL  JET,  by  J.B.Keller. 

Geer  and  Strikwerda  devised  a  slender  jet  theory  for  a  jet  falling  vertically,  with 
surface  tension  present.  They  solved  the  resulting  problem  numerically  by  extending  a 
method  they  had  developed  previously.  Their  results  showed  that  the  cross-section  of  the 
jet  decreased  in  area  and  oscillated  in  shape  with  distance  along  the  jet.  They  compared 
the  oscillations  with  Rayleigh’s  results  for  oscillations  of  a  jet  of  constant  circular  cross- 
section.  There  was  agreement  qualitatively  but  not  quantitatively. 

We  shall  analyze  the  oscillations  of  a  vertically  falling  jet  on  the  basis  of  their  theory 
in  order  to  obtain  a  quantitative  description  of  them  which  agrees  with  the  numerical 


results.  First  we  find  the  solution  for  a  vertical  jet  with  a  circular  cross-section.  Then  we 
determine  its  small  amplitude  oscillations.  Our  analysis  dillcrs  from  Rayleigh’s  because 
our  unperturbed  jet  is  falling,  so  its  velocity  aud  radius  arc  not  constant. 


4.  SURVIVAL  ESTIMATION  WITH  CENSORED  DATA, by  A.  Whittemore  and  J.B.  Keller 

New  nonparametric  methods  are  given  for  estimating  survival  probability  using  ran¬ 
domly  right- censored  data.  A  class  of  estimators  is  obtained  by  the  maximum  likelihood 
method,  in  which  the  hazard  rate  is  approximated  by  a  suitably  chosen  spline  function. 
The  class  includes  the  estimator  proposed  by  Nelson  (1969)  and  Altshuler  (1970).  The 
estimators  are  shown  to  be  uniformly  consistent  and  to  have  the  same  asymptotic  weak 
convergence  properties  as  the  Kaplan- Meier  estimator.  In  small  and  in  heavily  censored 
samples,  new  estimators  in  the  class  have  uniformly  smaller  mean  squared  error  than  do  the 
Kaplan- Meier  (1958)  and  Nelson-Altshuler  estimators,  according  to  computer  s;  lations. 
The  methods  are  extended  to  provide  new  estimates  for  both  the  baseline  ard  rate 
and  the  regression  coefficients  in  the  proportional  hazards  model  proposed  by  .*  (1972). 
Several  existing  estimates  for  these  quantities,  including  that  obtained  using  partial 
likelihood,  occur  as  special  cases  of  the  general  procedure.  The  methods  ar>  ..\  rated 
using  experimental  carcinogenesis  data. 


5.  EFFECTIVE  VISCOSITY  OF  A  PERIODIC  SUSPENSION,  by  K.Nunan  and  J.B.  Keller 

The  effective  viscosity  of  a  suspension  is  defined  to  be  the  four-tensor  which  relates  the 
average  deviatoric  stress  to  the  average  rate  of  strain.  We  determine  the  effective  viscosity 
of  an  array  of  spheres  centered  on  the  points  of  a  periodic  lattice  in  an  incompressible 
Newtonian  fluid.  The  formulation  involves  the  traction  exerted  on  a  single  sphere  by 
the  fluid,  and  an  integral  equation  for  this  traction  is  derived.  For  lattices  with  cubic 
symmetry  the  effective  viscosity  tensor  involves  just  two  parameters.  They  are  computed 
numerically  for  simple,  body-centered  and  face-centered  cubic  lattices  of  spheres  with 
solute  concentrations  up  to  90%  of  the  close-packing  concentration.  Asymptotic  results 
for  high  concentrations  are  obtained  for  arbitrary  lattice  geometries,  and  found  to  be 
in  good  agreement  with  the  numerical  results  for  cubic  lattices.  The  low  concentration 
asymptotic  expansions  of  Zuzovsky  also  agree  well  with  the  numerical  results. 
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6.  THE  ZERO  DISPERSION  LIMIT  OF  THE  KORTEWEG-DE  VRIES  EQUATION 
FOR  INITIAL  POTENTIALS  WITH  NON-TRIVIA L  REFLECTION  COEFFR  TENT, 
by  S.  Ven  abides. 


The  inverse  scattering  method  is  used  to  determine  the  distribution  limit  as  i  —  >  0  of 
the  solution  u(x,t,  e)  of  the  initial  value  problem: 

ut  —  6uux  +  e2uxxx  =  0 

u( x,  0 )  =  v(x) 

where  v(x)  is  a  positive  bump  which  decays  sufficiently  fast  as  x  — *  ±cc.  The  case  v[x)  <  0 
has  been  solved  by  Peter  D.  Lax  and  C.D.  Levermore. 

The  computation  of  the  distribution  limit  of  it(x,t,e)  as  e  —*  0  is  reduced  to  a 

quadratic  maximization  problem. 

7.  VELOCITY  AND  DENSITY  OF  A  TWO-DIMENSIONAL  ACOUSTIC  MEDIUM 
FROM  POINT  SOURCE  SURFACE  DATA,  by  S.  Coen,  Margaret  Cheney  and  A.  Weglein. 

An  inverse  acoustic  scattering  theory  and  algorithm  is  presented  for  the  reconstruction 
of  a  two  dimensional  inhomogeneous  acoustic  medium  from  surface  measurements.  The 
measurements  of  the  surface  pressure  due  to  a  harmonically  oscillating  surface  point  source 
at  two  arbitrary  frequencies  allows  the  separate  reconstruction  of  the  density  and  velocity 
of  the  subsurface.  This  is  a  first  step  towards  solving  the  inverse  problem  of  exploration 
geophysics. 


8.  GENETIC  VARIABILITY  DUE  TO  GEOGRAPHICAL  INH OMOGENEi f  Y,  by  J.B. 
Keller. 

The  frequency  ratio  of  two  alleles  is  studied  as  a  function  of  position  and  time,  in  a  two- 
dimensional  region,  by  means  of  a  nonlinear  diffusion  equation.  Each  allele  is  assumed  to 
have  a  selective  advantage  in  some  part  of  the  region.  An  asymptotic  solution  is  constructed 
when  the  selection  coefficient  is  large  compared  to  the  diffusion  coefficient,  be.,  when 
selection  acts  more  rapidly  than  diffusion.  As  time  increases,  the  solution  tends  to  an 
equilibrium  distribution  in  which  both  alleles  are  present  everywhere,  each  predominating 
where  it  has  the  advantage. 


